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The effective lifetime of the meson in the decay mode — t- J/^/o(980) is 
measured using 1.0 fb~^ of data collected in pp collisions at ^/s = 7TeV with the 
LHCb detector. The result is 1.700 ± 0.040 ± 0.026 ps where the first uncertainty 

7—i is statistical and the second systematic. As the final state is CP-odd, and CP 

^ violation in this mode is measured to be small, the lifetime measurement can be 

translated into a measurement of the decay width of the heavy B^ mass eigenstate, 

%J Th = 0.588 ±0.014 ±0.009 ps-i. 
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The decay J/^/o(980), /o(980) ^ tt+tt-, discovered by LHCb [Ij at close to the 

predicted rate j2j, is important for CP violation js] and lifetime studies. In this Letter, we 
make a precise determination of the lifetime. The J/ip fo{980) final state is CP-odd, and 
in the absence of CP violation, can be produced only by the decay of the heavy (H), and 
not by the light (L), mass eigenstate j^. As the measured CP violation in this final 
state is small [s], a measurement of the effective lifetime, Tj/^f^, can be translated into a 
measurement of the decay width, Fh. This helps to determine the decay width difference, 
ATs = Fl — Fh, a number of considerable interest for studies of physics beyond the 
Standard Model [6]. Furthermore, this measurement can be used as a constraint in the fit 
that determines the mixing-induced CP- violating phase in decays, 0s, using the J/ijj4> 
and J/-ip fo{980) final states, and thus improve the accuracy of the (f)s determination [5|[7j. 

In the Standard Model (SM) 0s = — 2arg J" J* , if sub-leading penguin contributions are 

L cs ^cb J 

neglected, where the Vij are the Cabibbo-Kobayashi-Maskawa matrix elements, which has 
a value of —0.036 [Jig Qgj^g rad [sj. Note that the LHCb measurement of 0^ [s] corresponds 
to a limit on cos0s greater than 0.99 at 95% confidence level, consistent with the SM 
prediction. 

The decay time evolution for the sum of B^ and B^ decays, via the b — )■ ccs tree 
amplitude, to a CP-odd final state, /_, is given by [9] 



F (5° ^ /-) + F (5° ^ /_) = ye-^=* <j e^r=*/2(l + cos0s) + e'^^^'/^^ 



cos 




where A/" is a time-independent normalisation factor and Fg is the average decay width. 
We measure the effective lifetime by describing the decay time distribution with a single 
exponential function 

F (5° ^ /_) + F (5° ^ /_) = Are"*/-^/^'/o . (2) 

Our procedure involves measuring the lifetime with respect to the well measured B^ 
lifetime, in the decay mode B^ ^ J/ipK*^, K*^ ^ K^ti^ (the inclusion of charge conjugate 
modes is implied throughout this Letter). In this ratio, the systematic uncertainties largely 
cancel. 

The data sample consists of L0fb~^ of integrated luminosity collected with the LHCb 



detector 10 in pp collisions at the LHC with 7 TeV centre-of-mass energy. The detector is 
a single- arm forward spectrometer covering the pseudorapidity range 2 < 77 < 5, designed 
for the study of particles containing 6 or c quarks. The detector includes a high precision 
tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction 
region, a large-area silicon-strip detector located upstream of a dipole magnet and three 
stations of silicon-strip detectors and straw drift-tubes placed downstream. Charged 
hadrons are identified using two ring-imaging Cherenkov detectors. Muons are identified 
by a muon system composed of alternating layers of iron and multiwire proportional 
chambers. The trigger consists of a hardware stage, based on information from the 
calorimeter and muon systems, followed by a software stage that applies a full event 
reconstruction. 
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Figure 1: Invariant mass distributions of selected (a) 7r"'"7r~ and (b) K~7r~^ combinations 
(solid histograms) for events within ±20 MeV of the respective and mass peaks. 
Backgrounds (dashed histograms) are determined by fitting the J/ipn^ir' {J/tp K~7i'~^) 
mass in bins of tt'^tt^ (K^tt^) mass. Regions between the arrows are used in the subsequent 
analysis. 



The selection criteria we use 
in -5°—)- J/ip7i^7i~ decays 



br this analysis are the same as those used to measure 
Events are triggered by a J/ijj — )■ decay and a 



11 

Boosted Decision Tree (BDT)~][l2] is used to set the J/%I)'k~^t^~ selection requirements 
As the effective i?° — t- J/?/'/o(980) lifetime is measured relative to that of the decay 
5° —7- J/ipK*^, we use the same trigger and BDT to select J/iIjK~tt~^ events, except for 
the hadron identification that is applied independently of the BDT. The selected tt+tt^ 
and K~n~^ invariant mass distributions, for candidates with J/t/'tt+tt" {J/ip K~7t'^) mass 
within ±20 MeV (we work in units where c = h = 1) oi the respective B mass peaks 
are shown in Fig. [T| The background distributions shown are determined by fitting the 
J/ipn^TT^ {J/ip 7c^) mass distribution in bins of tt+tt^ {K^ti^} mass. Further selections 
of ±90 MeV around the /o(980) mass and ±100 MeV around the K*^ mass are applied. 
The /o(980) selection results in a 5° — J/ip fo{980) sample that is greater than 99.4% 
CP-odd at 95% confidence level |13]. 

The analysis exploits the fact that the kinematic properties of the B^—^ J/^/'/o(980) 
decay are very similar to those of the B^ — t- J/ip K*^ decay. We can select B mesons in 
either channel using identical kinematic constraints and hence the decay time acceptance 
introduced by the trigger, reconstruction and selection requirements should almost cancel in 
the ratio of the decay time distributions. Therefore, we can determine the 5° — )■ J/ip /o(980) 
lifetime, tj/^j^^, relative to the — )• J/ifjK*^ lifetime, Tj/^^*o, from the variation of the 
ratio of the B meson yields with decay time 

R{t) = i?(0)e"*^^/"'^/"^-fo-i/"^/'^^*°) = i?(0)e-*^-^/*/o , (3) 
where the width difference Aj/^/^ = 1/rj/^/o - ^/tj/^k*o- 
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tips) 



Figure 2: Ratio of decay time acceptances obtained from simulation between — )■ 
J/ip fo{980) and 5°— >■ J/ipK*^ decays. The solid (blue) line shows the result of a linear fit. 



We test the cancellation of acceptance effects using simulated -8°—)- 7/^/^/0(980) and 
— 7- J/ipK*^ events generated using Pythia 6.4 |14 with a specific LHCb configura- 



tion (15], where decays of hadronic particles are described by EvtGen [Tg], and the LHCb 



detector simulation 17 based on Geant4 18 . Due to the selection requirements, the 



acceptances for both decay modes are equal to at t = 0, after which there is a sharp 
increase, followed by a slow variation for t greater then 1 ps. Based on this, we only use 
events with t greater than 1 ps in the analysis. To good approximation, the acceptance 
ratio is linear between 1 and 7ps, with a slope of a = 0.0125 ± 0.0036 ps~^ (see Fig. [2 
We use this slope as a correction to Eq. [3] when fitting the measured decay time ratio 

R{t) = i?o(l + at)e-^^-'/^fo. (4) 

Differences between the decay time resolutions of the decay modes could affect the 
decay time ratio. To measure the decay time resolution, we use prompt events containing 
a J/ip meson. Such events are found using a dimuon trigger without track impact 
parameter requirements, plus two opposite-charged tracks with similar selection criteria 
as for J/iIjtt~^7c~ (J/ip K~7i~^) events, apart from track impact parameter requirements, 
including that the J/tp7r^7T~ [J/ip K~7t^) mass be within ±20 MeV of the (B^) mass. 
These events are dominated by zero lifetime background. We use a triple-Gaussian 
function for the resolution model and find the average effective decay time resolution for 
J/V^/o(980) and 5°^ J/'^K*'^ decays to be 41.0±0.9fs and 44.1±0.2fs respectively, 
where the uncertainties are statistical only. This difference was found not to bias the decay 
time ratio using simulated experiments. 

In order to determine the — )■ J/ip fo{980) lifetime, we determine the yield of B 
mesons for both decay modes using unbinned maximum likelihood fits to the B mass 
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Figure 3: Invariant mass distributions of selected (a) J/^/'tt+tt" and (b) J/iI)K~tt~ candi- 
dates. The solid (blue) curves show the total fits, the long dashed (purple) curves show 
the respective -8°— )■ 7/^/^/0(980) and )■ J/ipK*^ signals, and the dotted (grey) curve 
shows the combinatorial background. In (a) the short dashed (light blue) curve shows the 
5°— )• J/ipn^n^ background and the dash dotted (green) curve shows the 5°—)- J/ipK^ tt"*" 
reflection. In (b) the short dashed (pink) curve shows the -Bg— )■ J/ipK^Ti^ background. 



distributions in 15 bins of decay time, of equal width between 1 and 7ps. We perform a 
fit to the ratio of the yields as a function of decay time and determine the relative lifetime 
according to Eq. |4| We obtain the signal and peaking background shape parameters 
by fitting the time-integrated dataset. In each decay time bin, we use these shapes and 
determine the combinatorial background parameters from the upper mass sidebands, 
5450 < m{J/iljfo) < 5600 MeV and 5450 < m{J/^K*^) < 5550 MeV. With this approach, 
the combinatorial backgrounds are re-evaluated in each bin and we make no assumptions 
on the shape of the background decay time distributions. This method was tested with 
high statistics simulated experiments and found to be unbiased. 

The time-integrated fits to the J/^/o(980) and the J/ipK*^ mass spectra are shown 
in Fig. |3] The signal distributions are described by double Crystal Ball functions (sums 
of two Crystal Ball functions \T9\ with equal means but tails on opposite sides) and 
the combinatorial backgrounds by exponential functions. All parameters are determined 
from data. There are 4040 ± 75 5° ^ J/V^/o(980) and 131920 ± 400 5° J/^K*^ 
signal decays. The decay time distributions, determined using fits to the invariant mass 
distributions in bins of decay time as described above, are shown in Fig. |4j These are 
made by placing the fitted signal yields at the average — )■ J/ip K*^ decay time within 
the bin rather than at the centre of the decay time bin. This procedure corrects for the 
exponential decrease of the decay time distributions across the bin. The subsequent decay 
time ratio distribution is shown in Fig. |5} and the fitted reciprocal lifetime difference is 
'^j/ipfo = —0.070 ± 0.014ps~^, where the uncertainty is statistical only. Taking Tj^^^^.o to 



be the mean 5° lifetime 1.519 ± 0.007 ps 20 , we determine tj/^j^ = 1.700 ± 0.040 ps 



4 



Oh 

o 



600 



400 



T 
ICQ 



200 





f (ps) 



f (ps) 



Figure 4: Decay time distributions for (a) )■ J/ip fo(980) and (b) 5°— > J/ipK*^. In (b) 
the error bars are smaller than the points. 




Figure 5: Decay time ratio between i?^ — )• J/-?/'/o(980) and — )• J/ipK*^, and the fit 
for Aj/^/,. 



Sources of systematic uncertainty on the -B°— )■ J/ip fo {980) lifetime are investigated and 
listed in Table [T| We first investigate our assumptions about the signal and combinatorial 
background mass shapes. The relative change of the determined J/ip foiQSO) lifetime 

between fits with double Crystal Ball functions and double Gaussian functions for the 
signal models is 0.001 ps, and between fits with exponential functions and straight lines 
for the combinatorial background models is 0.010 ps. The different particle identification 
criteria used to select 5° J/ijj fo{980) jj,^ fi^n^n' and i?° J/tjjK*^ -)■ /x+/x~i^'~7r+ 
decays could affect the acceptance cancellation between the modes. In order to investigate 
this effect, we loosen and tighten the paxticle identification selection for the kaon and repeat 
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Table 1: Summary of systematic uncertainties on the -B^— )■ J/ip fo{980) effective lifetime. 



Source Uncertainty (ps) 

Signal mass shape 0.001 

Background mass shape 0.010 

Kaon identification 0.007 

Acceptance 0.018 

Statistical bias 0.012 

CP-even component 0.001 



lifetime 20 



0.009 



Sum in quadrature 0.026 



the analysis. The larger difference with respect to the default selection, 0.007 ps, is assigned 
as a systematic uncertainty. We also assign half of the relative change between the fit 
without the acceptance correction and the default fit, 0.018 ps, as a systematic uncertainty. 
Potential statistical biases of our method were evaluated with simulated experiments using 
similar sample sizes to those in data. An average bias of 0.012 ps is seen and included 
as a systematic uncertainty. The observed bias vanishes in simulated experiments with 
large sample sizes. As a cross-check, the analysis is performed with various decay time 
bins widths and fit ranges, and consistent results are obtained. The possible CP-even 



component, limited to be less than 0.6% at 95% confidence level 13 , introduces a 0.001 ps 
systematic uncertainty. Using the PDG value for the lifetime as input requires the 
propagation of its error as a systematic uncertainty. All the contributions are added in 
quadrature and yield a total systematic uncertainty on the lifetime of 0.026 ps (1.5%). 
Thus the effective lifetime of the J/ip fo(980) final state in 5° decays, when describing the 
decay time distribution as a single exponential is 

Tj/^^^ = 1.700 ± 0.040 ± 0.026 ps . (5) 

Given that (f)s is measured to be small, and the decay is given by a pure b — )■ ccs tree 
amplitude, we may interpret the inverse of the -B^ — )■ J/ip fo{980) effective lifetime as a 
measurement of Fh with an additional source of systematic uncertainty due to a possible 
non-zero value of 0s. For cos 0s = 0.99, Fg = 0.6580 ps~^ and ATg = 0.116 ps~^ [5], tj/^j^ 
changes by 0.002 ps. This is added in quadrature to the systematic uncertainties on Tj/^f^ 
to obtain the final systematic uncertainty on Fh. 

In summary, the effective lifetime of the 5° meson in the CP-odd J/^/o(980) final 
state has been measured with respect to the well measured lifetime in the final state 
J/ipK*^. The analysis exploits the kinematic similarities between the -Bg— )■ J/ip fo{980) 
and B^ — )■ J/tp K*^ decays to determine an effective lifetime of 

Tj/^/o = 1.700 ± 0.040 ± 0.026 ps. 
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corresponding to a width difference of 

Aj/^/o = -0.070 ±0.014 ±0.001 ps"\ 

where the uncertainties are statistical and systematic respectively. This result is consistent 
with, and more precise than, the previous measurement of 1.70 1*^0 ±0.03 ps from CDF 
Interpreting this as the lifetime of the heavy eigenstate, we obtain 

Th = 0.588 ± 0.014 ± 0.009 ps~\ 

This value of Fh is consistent with the value 0.600 ± 0.013 ps~^, calculated from the values 
of Fo and AF., in Ref. |5l . 
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